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The Helfand-Werthamer (HW) scheme^ of evaluating the orbital upper critical field is generalized 
to anisotropic superconductors in general, and to two-band clean materials, in particular. Our formal 
procedure diflfers from those in the literature; it reproduces not only the isotropic HW limit, but also 
the results of calculations for the two-band superconducting MgB2'^'^ along with the existing data on 
Hc2{T) and its anisotropy "f{T) — Hc2,ab{T) / Hc2,c{T) (a, c are the principal directions of a uniaxial 
crystal). Using rotational ellipsoids as model Fermi surfaces we apply the formalism developed to 
study 7(r) for a few different anisotropies of the Fermi surface and of the order parameters. We find 
that even for a single band d-wave order parameter 'y(T) decreases on warming, however, relatively 
weakly. For order parameters of the form A{kz) = Ao(l + ri cos kztt),* according to our simulations 
7(r) may either increase or decrease on warming even for a single band depending on the sign of 
T]. Hence, the common belief that the multi-band Fermi surface is responsible for the temperature 
variation of 7 is proven incorrect. 

For two s-wave gaps, 7 decreases on warming for all Fermi shapes examined. For two order 
parameters of the form A{kz) = Ao{l + 7} cos kza) presumably relevant for pnictides, we obtain •y{T) 
increasing on warming provided both r?i and r?2 are negative, whereas for rj's > 0, 7(T) decreases. 
We study the ratio of the two order parameters at Hc2{T) and find that the ratio of the small gap 
to the large one does not vanish at any temperature even at Hc2{T), an indication that this does 
not happen at lower fields. 



I. INTRODUCTION 

The seminal work of Helfand and Werthamer (HW)^ 

on the temperature dependence of the upper critical field 
Hc2{T) is routinely applied to analyze data on new su- 
perconductors despite the fact that HW considered the 
isotropic s-wave case whereas in themajority of new ma- 
terials both Fermi surfaces and order parameters are 
strongly anisotropic. The problem of Hc2{T) has been 
studied theoretically for anisotropic situations as well: for 
layered systems,^ for a few cases of hexagonal anisotropy 
of the Fermi surface and of the order parameter,® for 
the two-band MgB2,^'^'''"^ for d-wave cuprates, see. e.g., 
Ref. 10 and references therein, and in a comprehensive 
work in Ref. 11, to name a few. 

A ubiquitous feature of Hc2 in many new superconduc- 
tors is the temperature dependent anisotropy parameter 
7 = Hc2.ab/ Hc2,c (for Uniaxial materials), the property 
absent in conventional one-band isotropic s-wave materi- 
als. For example, 7(T) of MgB2 decreases on warming, 
whereas for many Fe-based materials j(T) increases with 
increasing T.^^^^^ Up to date, the T dependence of the 
anisotropy parameter 7 is considered by many to be 
caused by a multi-band character of the materials in ques- 
tion with the commonly given reference to the example 
of MgB2. To our knowledge, no explanation was offered 
for the "unusual" j{T) of pnictides. 

In this work we develop a method to evaluate both 
Hc2,c{T) and 7(T) that can be applied with minor mod- 
ifications to various situations of different order parame- 
ter symmetries and Fermi surfaces, two bands included. 

Having in mind possible applications for data analysis, 
we provide only the necessary minimum of analytic de- 



velopments resorting to numerical methods as soon as 
possible. 

The upper critical field is affected by many factors: 
magnetic structures that may coexist or interfere with su- 
perconductivity, the paramagnetic limit, scattering, etc. 
In this work we have in mind to establish a qualita- 
tive picture of how the general features of anisotropic 
Fermi surfaces and order parameters affect i?c2 and its 
anisotropy, in particular. 

Since the paramagnetic effects and the possibility of 
Fulde-Ferrel-Larkin-Ovchinnikov phase are not included 

in our scheme (one can find a comprehensive discussion of 
these questions in Ref. 8), applications to materials with 
high Hc2{0) should be carried out with care; our results 
can prove useful for interpretation of data at tempera- 
tures where H^2{T) does not exceed the paramagnetic 
limit. 

As far as applications to two-band materials are con- 
cerned, we note that our formalism applies only to super- 
conductors with a single critical temperature Tg. More 
exotic possibilities of two component condensates with 
two distinct T^'s are out of the scope of this paper; these 
are considered on general group-theoretical symmetry 
grounds, e.g., in Refs. 16 and 17. 

Fine details of Fermi surfaces are unlikely to strongly 
influence Hc2(T) because, as is well-known and shown 
explicitly below, only the integrals over the whole Fermi 
surface enter equations for Hc2{T). Circumstantial evi- 
dence for a weak connection between fine particularities 
of the Fermi surface and Hc2{T) is provided by the very 
fact that the HW isotropic model works so well for many 
one-band s-wave materials, although their Fermi surfaces 
hardly resemble a sphere, take e.g. Nb. Another exam- 



2 



pie is given by the calculations of Refs. 3 and 2 for MgB2 
based on different model Fermi surfaces, but giving simi- 
lar results reasonably close to the measured Hc2{T). We, 
therefore, model actual Fermi surfaces of uniaxial mate- 
rials of interest here by rotational ellipsoids (spheroids) 
choosing them as to have averaged squared Fermi veloc- 
ities equal to the measured values or to those calculated 
for realistic Fermi surfaces. This idea, in fact, has been 
employed by Miranovic et al. for MgB2.^ Also, we tested 
our method on a rotational hyperboloid as an example 
of open Fermi surfaces, Appendix D. This work is still in 
progress. We intend to study more about effects of open 
Fermi surfaces (or two-band combinations of closed and 
open surfaces) on Hc2{T). 

We focus in this work on the clean limit for two major 
reasons. First, commonly after discovery of a new su- 
perconductor, an effort is made to obtain as clean single 
crystals as possible since those provide a better chance to 
study the underlying physics. Second, almost as a rule, 
new materials are midti-band so that the characteriza- 
tion of scattering leads to a multitude of scattering pa- 
rameters which cannot be easily controlled or separately 
measured. Besides, in general, scattering suppresses the 
anisotropy of ffc2 , the central quantity of interest in this 
work. We refer readers interested in effects of scattering 
to a number of successful dirty-limit microscopic models, 
e.g. to Ref. 7 and references therein. 

In the next Section, we begin with the general dis- 
cussion of the Hc2 problem for arbitrary Fermi surfaces 
and order parameters. We show that in the isotopic 
limit our approach is reduced to that of HW. The basic 
HW approach is then applied to anisotropic situations. 
The derivation involves rescaling of the coordinates and, 
therefore, necessitates co- and contravariant vectors rep- 
resentations. In our view, disregarding this necessity may 
lead to incorrect conclusions. Another formal feature of 
our approach should be mentioned: wc avoid the mini- 
mization relative to the actual coordinate dependence of 
the order parameter in the mixed state often employed 
to find Hc2{T)?'^''^'^^'^^ In this sense, our method is close 
to the original HW work that stresses that Hc2 {T) is ac- 
tually an eigenvalue problem. 

Next, we formulate the problem for two s-wave bands 
and show that along with Hc2 (T) and its anisotropy one 
can find the ratio of order parameters on two bands, a 
quantity that up to date has been studied only in zero- 
field. We then formulate the problem for two bands with 
order parametcirs of different symmetries. To show how 
the method actually works, we consider in detail one or 
two bands with Fermi surfaces as rotational ellipsoids. 
The method is demonstrated on the well-studied example 
of MgB2. 

The anisotropy parameter 7 is shown to depend on 
temperature even for the one-band case for other than 
s-wave order parameters. This dependence is weak the 
d-wave materials with closed Fermi surfaces, is stronger 
for open Fermi shapes, and is stronger yet for order pa- 
rameters of the form A = Ao(l -|- r/cosfc^a), one of the 



candidates for Fe-bascd materials.^ Moreover, 7(T) in- 
creases or decreases on warming depending on the sign 
of the coefficient r], in other words, on whether A is max- 
imum or minimum at kz = 0. These results challenge 
the common belief that temperature dependence of 7 is 
always related to the multi-band topology of Fermi sur- 
faces. 

For two bands, after checking the method on MgB2, 
we focus on situations with dominant inter-band cou- 
pling which is relevant for Fe-based materials. While in 
most cases we have considered, Hc2,c{T) is qualitatively 
similar to the HW curve, the anisotropy parameter 7(T) 
may show a non-monotonic T dependence even for s-wave 
order parameters depending on the Fermi surface shapes 
and densities of states (DOS). Most interesting are the or- 
der parameters A = Ao(l + 77 cos kzu) which yield nearly 
linear increase of j{T) on warming, a ubiquitous feature 
seen in many Fe-based superconductors. 

II. THE PROBLEM OF Ha2 

Our approach is basically that of HW, although for- 
mally the equations we solve for Hc2 {T) are different and 
can be applied to anisotropic and multi-band situations. 

To establish the link to HW, we start the discussion with 
the one-band case. The problem of the 2nd order phase 
transition at Hc2 is addressed here on the basis of lin- 
earized (the order parameter A ^ 0) quasiclassic Eilen- 
berger equations. For clean materials we have: 

{2L0 + v-U)f = 2A/h, (1) 

Here, v is the Fermi velocity, 11 = V + 2TTiA/(pQ, A is the 
vector potential, and (po is the flux quantum. A(r,feir) 
is the order parameter that in general depends on the 
position kp at the Fermi surface (or on v). The function 
fir.v.Lu) originates from Gor'kov's Green's function in- 
tegrated over energy near the Fermi surface. Further, 
A''(0) is the total density of states at the Fermi level per 
spin; the Matsubara frequencies are lo = TTT{2n + l)/h 
with an integer n. The averages over the Fermi sur- 
face are shown as (...). The Eilenberger function g = 
\/l — ff^ = 1 at Hc2- The temperature T is in energy 
units, i.e., fcs = 1. 

The self-consistency equation (2) is written for the gen- 
eral case of anisotropic gaps: A = ^'(r, T) $l(kir). The 
function ri(ki?) which determines the kp dependence of 
A is normalized so that 

{n') = i, (3) 

for details see, e.g., Ref. 20. Eq. (2) corresponds to the 
factorizable coupling potential, V{k,k') = VoO,{k)Cl{k'). 

This popular approximation^-'^ works well for one band 
materials with anisotropic coupling. We show in Sections 
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V and VI how this convenient shortcut can be generalized 
to a multi-band case. 

We now recast the self-consistency Eq. (2) by writing 
the solution of Eq. (1) in the form 



with 



9 r°° 
" Jo 

by using the identity 



1 



2 f°° 



and by summing up over to: 

Jo sinh u\ / 



(4) 



(5) 



(6) 



where u = 2TTTp/h and t = T/Tc- Hence, wc got rid of 
the summation over w, a convenient feature for further 
analysis. 

The self-consistency Eq. (6) can be further rewritten in 
the form free of the divergent factor 1 / sinh u in the inte- 
grand. To this end, we integrate by parts the right-hand 
side (RHS) of Eq. (6) using du/sinhu = d lntanh(u/2). 
The first term on the RHS diverges: 



f 

Jo 



du 



-^Intanh 



TrTp 



/o sinh u 
The second term transforms to: 



p->-0 



(7) 



Intanh^e-^'^n 



p=0 



irTp 



rfplntanh^e"''^"*), 



(8) 



The first term here and the contribution (7) cancel, and 
we obtain instead of Eq. (6): 

*lni= / dp liit&nh^ {n'^vUe-P''-^^) . (9) 
Jo " 

Here, the singularity at p ^ is integrable. 

A. Hc2 near Tc 

In this domain, the gradients 11 ~ — >• 0, and one 
can expand exp(— pwll) in the integrand (9) and keep 
only the linear term: 

-^6t=^,^{ifivnr^), (10) 

where dxxlntanhx = -7C(3)/16 with C(3) = 1.202. 
This is, in fact, the anisotropic version of the linearized 
Ginzburg-Landau (GL) equation 



(11) 



{Vl?ViVk) , 



T=l-t, 



(12) 



the result of Gor'kov and Melik-Barkhudarov.^^ Solving 
the eigenvalue problem for Eq. (11) which is similar to 
one for a charged particle in uniform magnetic field, see 
e.g. work by Tilley,^^ one finds the critical fields in two 
principal directions of uniaxial materials: 



Hc2,c = 



so that 



H. 



c2,a 



H, 



c2,c 



Sag 
Sec 



The angular dependence 



Hc2{0) 



H. 



c2,a 



\J sin^ ^ + 7^ cos^ Q 



(13) 



(14) 



(15) 



is a direct consequence of Eq. (11) in which £|j, is a sec- 
ond rank tensor {d is the angle between the applied 
field and the c axis). We argue below that Eq. (11) 
holds at iJc2, in fact, at all temperatures, and so should 
the angular dependence (15), the common practice to 
call it "GL" notwithstanding. We show that this an- 
gular dependence holds for any order parameter sym- 
metry and for any Fermi surface shape including two- 
band situations. These conclusions are, in fact, confirmed 
experimentally^"^"^" and by calculations of Ref. 2. 



III. ISOTROPIC GAP ON A FERMI SPHERE 

This problem has been solved by HW for the whole 
curve Hc2{T)} It is instructive and useful for the follow- 
ing generalization to the anisotropic case to reproduce 
their result within the quasiclassic scheme. ■^^ 

It was established in Ref. 1 that at Hc2 {T) at any tem- 
perature, the order parameter satisfies a linear equation 



(16) 



in which £(T) should be determined so as to satisfy the 
self-consistency equation. One can see that this equa- 
tion is equivalent to the Schroedinger equation for a 
charge moving in uniform magnetic field and that the 
max;imum field in which non-trivial solutions A exist is 
Hc2 = 0o/27r£2. For the field along z we choose the gauge 
Ay = Hx. One readily verifies that in terms of operators 

U"^ =Ilx±iUy, Ua; = dx, Uy = dy + iq'^x, 

= 2nH/cPo , (17) 
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Eq. (16) reads H+R-A = provided = There- 
fore, we obtain a useful property of A at Hc2{T): 



where C = 0.577 is the Euler constant. The averaging 
over the Fermi sphere with vj_ — vp sin 9 is readily done: 



n-A = n-* = 0. 



(18) 



„C-2 



27r2T2 



We now introduce = v^iivy so that vU = {v 11+ - 



w+n )/2 and evaluate the average (^vlle p""^ ^! 

needed in the self-consistency Eq. (9). To this end, we 
use the known property of exponential operators: 



Thus, we obtain: 



27rC2(0) 



,2-C 



(26) 



(27) 



(19) 



Here, P = —pv n+/2, Q = —pv~^Il /2, the commutator 
[Q, P]/2 = -p^vl/4e, and vl=vl + v^. 
Since n~* = and e'^* = we have: 



2 2 



i-pr-u+y 



2"n! 



the value obtained variationally by Gor'kov'^^ and proven 
to be exact by HW. 

Close to Tc, Hc2 is obtained as isotropic limits of 
Eqs. (13). It is instructive, however, to see how this 
can be deduced directly from Eq. (24). In this domain, 
Int « t — 1 = —T and h oc t. Then, the integral in 
Eq. (24) should be evaluated in zero order in r, in other 
words, h in ciqy{— phs'^) can be set zero. Integration over 
s gives -(7C(3)/16)(Ai), whereas (p) = {vj_/vl) = 2/3: 



(20) 



/i = 12r/7C(3). 



(28) 



After averaging over the Fermi sphere, only {v~^v ) sur- The reduced HW field 



vives (use = v±e 



(21) 



h*{0) = 



i?c2(0) 



h{0) 7C(3)e2 



0.727, 



(29) 



Now, with the help of Eq. (16), the self-consistency 
Eq. (9) yields: 



^2 hit 



dp p In tanh 



irTp 



1 



2./o h 



Introducing the dimensionless field 



h = 



47r2T2e 



2 ■ 



(e-^vl) . (22) 



(23) 



and a variable s — nTcp/h, we rewrite Eq. (22) for ^{t) 
as an equation for h{t): 

POO . . 

lnt = 2hj dss Intanh(st) (^pe-'^'''J , (24) 

where IJ, = v'^/vp and Vp is the Fermi velocity. The field 
h is the upper critical field in units of (j)^/ {h?v'^p /2ttT'^). 
Thus, solving Eq. (24) with respect to h{t) we have the 
HW solution of the problem for the isotropic one-band 
case. In the following we refer to Eq. (24) as the HW 
result although in the original work they obtained a dif- 
ferent form of the equation for Hc2{T). 

At arbitrary T, Eq. (24) can be solved numerically; the 
exceptions are T — > and T — > Tc. Since /i(0) is finite, 
the integral over s is truncated at s ~ l/y^/xTi; therefore, 
for low enough t we have In tanh(st) w In f -|- In s. The 
integration over s then yields: 



\K2iTc)\Tc \h'{l)\ 48eC 
the correct HW value. 

IV. ANISOTROPIC ONE-BAND CASE 



The central point of the HW paper is the proof that the 
linearized GL equation (16) holds not only near but 
along the whole curve Hc2 (T) . For anisotropic materials 
Eq. (16) should be replaced with Eq. (11) where all com- 
ponents of the tensor {£,^)ik should be determined from 
the self-consistency equation. We consider here uniaxial 
materials for which the symmetric tensor has two in- 
dependent eigenvalues. One has for the field along one 
of the principal crystal directions which we call z: 



-(eLn^ + Cn')A = A. 



(30) 



We denote £,xx = ^/y/rrh,:, £.yy = H 
with dimensionless constants m 
quantities, ^xx, £,yy, and ^, 
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x.y,z- Since the three 
arc replaced with four, ^ and 



^x,y,z^ 



we can impose an extra condition: nixmymz 



1. For uniaxial materials of interest here with rric ^ 
nT-a = "i6i we introduce the anisotropy parameter 7 = 
sjmclnia^ so that all "masses" are expressed in terms 



of 



1 



-2/3 



C 

2" 



In 



27rr,e(0) 



(25) 



= 7 It is worth noting 
that rn's here do not necessarily have the meaning of 
the band theory effective masses; rather, they are pa- 
rameters describing the anisotropy of Hc2', near Tc they 
can be expressed in terms of Fermi velocities and the gap 
anisotropy, Eq(14).22 Hence, if the ansatz (11) is correct, 
the self-consistency equation must provide equations to 
determine ^(T) and j{T). 
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To make use of the property (18) in anisotropic situa- 
tions we rescale coordinates in Eq. (30): 

(31) 
(32) 

(33) 



d 



27tH 



dx' 2 ' \dy' 



2 A 



Formally, Eq. (32) is equivalent to the isotropic Eq. (16). 
The upper critical field is then determined by g'^ = 1/^^: 



c2 



277^2 • 

Therefore, we have in the uniaxial case: 



(34) 



H, 



c2,c 



277^2 



H. 



c2,b 



27r^' 



:\frrQmc (35) 



and 



c2.b 




7- 



(36) 



Thus, the formally introduced "masses" are related to 
the measurable ratio of Hc2^s. 

Any coordinate transformation results in transforma- 
tion of vectors (and tensors). The scaling transformation 
(31) necessitates the co- and contravariant representa- 
tions for vectors, see, e.g., Ref. 33 or 34. The covari- 
ant gradients tt^ = d/dx' and iTy = d/dy' + iq''^x' have 
the same properties as their isotropic counterparts ^x,v 
Eq. (32) acquires the "isotropic" form: 

-ei7rl + 7Tl)A = A, (37) 

TT-TT+A = -2A/^2^ 7r-A = 0. (38) 

The following manipulation then is similar to what was 

done in the isotropic case; one, however, should keep 
track of differences between co- and contravariant com- 
ponents of vectors. 

Since the scalar products are invariant, we have vTl = 
v^TTi = {iy~TT~^ + v'^'K~)/2, v'^ — ± iv^ where 

= Vx^m,x , = Vy^niy (39) 

are the contravariant components of the Fermi velocity 
that transform as coordinates, Eq. (31). Further, we use 
the property (19) of exponential operators with 



4^2 ■ 



and 



Since tt = and e'^^ = we have: 



-{v Tr'^ + v'^TT ) 



(-pz/-7r+)" 
2"n! 



n=0 



(40) 



(41) 



(42) 



where rj = p^v^/A^"^. 

The next step in the "isotropic derivation" was to use 
the fact that = 2?;|,/3 whereas all other averages 

(such as {v^v~v~)) vanish because v'^ = v±e^^'^ where 
(j) is the azimuthal angle on a sphere. To prove rigorously 
that this is valid for a general Fermi surface is difficult, 
although it is probably true for uniaxial materials of in- 
terest here. As mentioned in the Introduction, to make 
progress in evaluation of Hc2 and its anisotropy we resort 
to modeling actual Fermi surfaces with spheroids. The 
rescaling employed above, in fact, transforms spheroids 
to "spheres" in rcscaled coordinates so that we can still 
claim that the only surviving product after averaging of 
Eq. (42) is (i/+z/~) and we obtain: 

vHe"^^"*^ ^ -|(l22e-''jy+j.-)7r-7r+*. (43) 

The self-consistency Eq. (9) now yields with the help of 
Eq. (38): 

^'^\nt=\( dpplntanh^(02e-''z/2). (44) 

This equation (written for a certain field orientation) con- 
tains two unknown functions, ^(T) and ^(T). Therefore, 
one has to write it for two field orientations: (a) for H 
along c with ij = p^v^/A^"^ and 

vl = maiyl + «2) = 7-2/3 (^;2 + y^) , (45) 

and (b) for H along h with 

v'i = ruavl + rucvl = ^'^/^{vl + j^v^) . (46) 

In principle, by solving the system of these two equations 
one can determine both ^(T) and 7(r), thus proving the 
correctness of the ansatz (11). 

Since the Fermi velocity is not a constant at anisotropic 
Fermi surfaces, we normalize velocities on some value vq 
for which we choose 



(47) 



where Ep is the Fermi energy and A''(0) is the total den- 
sity of states at the Fermi level per spin. One easily 
verifies that vq = vp for the isotropic case. 

We now write the self-consistency Eq. (44) for H || c, 
i.e. with of Eq. (45), in dimensionless form. To this 
end, we go to the integration variable s = nTc/h, divide 
both parts by ^2^ and multiply and divide the integrand 
by v^: 

Inf = 2/ic / s Intanh(st) Ul^ Hce'^^-'^^'J ds , (48) 

, _ fi2^2^-2/3 



47r2T2^2 ' 



f 2 -I- 

^ 



(49) 



One can see that he is in fact Hc2.c in units of 
4>qI{T?v1I2-kTI). An important feature of Eq. (48) 
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should be noted: it docs not contain the anisotropy 7 
exphcitly so that it can be solved for hc{t). 

Writing Eq. (44) for H ± c with Uj_ of Eq. (46) we 
obtain: 

Ini = 2/ic / s Intanh(si) /o^be"'"'''"*^) ds , (50) 
= -o • (51) 



For the isotropic s-wavc case, 7 = f7 = 1 and Eqs. (48) 
and (50) coincide with each other and with Eq.(24) of 
HW. We note that the integrals on the RHS of Eqs. (48), 
(50) differ only in /x's; for brevity we denote: 



s In tanh(st) (o ^e""'*"*' ) ds . (52) 



1 = 



Thus, the general scheme for solving for Hc2{T) and its 
anisotropy consists of (a) solving Eq. (48) for hc{t) and 
(b) for the now known hcit), solving Eq. (50) for ^{t). 



A. T Tc and T 

Analysis of Eqs. (48), (50) for arbitrary temperatures 
is difficult because hc{t) and t enter the integrals X in 
a nonlinear fashion. The situation is simpler near 
where Int « — r = t — 1 and he cx t. Therefore, X can 
be evaluated in zero order in r, in other words, he in the 
exp{— fi^hcX^) can be set zero and t = 1 in lntanh(a;t). 
We obtain after integration: 



1 = 



16 \ 



r2> 



(53) 



where for iJ || c one takes fi = Hc whereas fi = fib for 
H ± c. Eq. (48) now yields: 

8t 

7C(3) {n^fic 

It is readily shown that Eq. (50) for j{Tc) reduces to 



hr = 



7C(3) {n^fic) 



(55) 



Using fib and fic of Eqs. (51), (49) one reproduces the 
general result (14). 

As t ^ 0, he ^ const, and the exponential factor in I 

truncates the integrand at a finite x ^ 1/^/Jih^. Hence, 
at small enough t, lntanh(a;<) ~ ln{xt) = Int + hnx. One 
can now integrate over x: 



X{i2) = 



Int C+{fl'^ln{hcfi)) 



2/ic 4/ic 
Substituting this in Eq. (48) one obtains: 

hc{0) = exp (-C - (f2^ In/ic)) • 
The HW ratio for a general anisotropy reads: 



MO) 



^(f2V)exp(-(02ln^,)) 



(56) 



(57) 



(58) 



Thus, the HW number h*{0) = 0.727 is corrected by 
both Fermi surface shape and by the order parameter 
symmetry. 

Eq. (50) for 7(0) along with expression for /ic(0) readily 
gives at T = 0: 



(n^lniie) = (n'^inub). 



TWO S-WAVE BANDS 



(59) 



The general self-consistency equation for two bands 
reads: 

A«(r,fe) = 27rT^iV^(y„/3(fc,fe')//3(r,fe',a;))^^. (60) 

0,u> 

Here, a,fi = 1,2 are band indices and are the bands 
densities of states (DOS'). We consider elements of Va^ 
as constants, so that our model is a weak-coupling two- 
band theory in which the s-wave (i.e., k independent) 
order parameters Ai.2(t', T, H) should be calculated self- 
consistently for a given coupling matrix Vq,^. 

As is commonly done, it is convenient to rewrite this 
equation in the form containing the measured critical 
temperature which is related to the effective coupling Vq 
via the standard BCS formula 



A(0) = 7re-^T, = 2;tt^^e-VJ^(o)^o ^ 



(61) 



where huju is the energy scale of a "glue boson". To 
this end, we introduce the normalized coupling matrix 
= Vafi/Vo and use the relation identical to Eq. (61) 
for T,: 



1 



NiO)Vo 
We then obtain 



In- 



T 



(62) 



w>0 



hhJ 



(63) 



with np = Np/N{G). 

Solving the self-consistency equation in zero field and 
A — one obtains a relation for Tg (or for Vq) in terms 
of couplings Va0 and DOS' Ua (Appendix B): 

V^'^nm^d - V'K^ii + ^2^22) + 1 = 0, (64) 
d = VxxVii - Vli , (65) 

Therefore, the normalized \aj) = Vajs/Vo obey 

nin2S — niAii — n2A22 -1- 1 = 0, (66) 
<5 = AiiA22-Af2. (67) 

This property, in fact, means that normalized couplings 

^ai3 for a given have only two independent compo- 
nents, which are chosen in the following as An and A22. 
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in literature: the latter are A^'* 

and A 5^*2 



To avoid misunderstanding, we stress that our notation 
for the normalized Ac,^ = V^/s/Vq differs from A^^ used 

NiVi2. It should also be noted that Eq. (63) 
is not valid for the unlikely situation of zero inter-band 
coupling, V12 = 0, because two decoupled condensates in 
general have two different critical temperatures. 

As was done in Section II, we can transform the self- 
consistency Eq. (63) to: 



and in Eq. (46) for H L c. Substituting this in the sys- 
tem (68) we obtain after straightforward algebra: 



aiiAi -I- ai2A2 = , 
a2iAi +a22A2 = 0, 



with 



A„lnt = ^nfiXa, 



■ [ dplntcaih^(vUe-P''^A) , ~ 2 L 

Jo n \ / /3 



On = ^ Int-AiiJi, ai2 = -A12 J2 , 
a2i = -A21 Ji , a22 = Int - A22 J2, 



. (j/oplntanh S lv\e 



(72) 

(73) 

(74) 



(68) 

the averaging in the last term here is only over /3-band. 

We have generalized the HW isotropic one-band 
approach by showing that linearized GL equation 
~(C^)i/£lIiIIfcA = A holds everywhere along the i?c2 line. 
Clearly, the tensor gives the length scales of spatial 

variations of A at H^^- Before solving the self-consistency 
equation for two-band systems (68) it is instructive to 
recall the situation in the GL domain of a two-band ma- 
terial, which has recently been discussed in some detail 
for the case of two isotropic bands. 

The system of two GL equations for two order param- 
eters written in terms of coefficients of the GL energy 
expansion looks-at first sight as containing two different 
coherence lengths, i.e., each order parameter varies in 
space with its own length scale different from the other. 
It has been shown, however, that at T = Tc these two 
length scales coincide, provided of course that the sys- 
tem has a single T^?"^ In fact, two GL equations can be 
written as one with a single coherence length ^ which is 
related in a non-trival manner to coefficients of the GL 
energy functional. Thus, for a material with two isotropic 
bands, the linearized GL equation is the same for both 
bands: 



-e^n^A„ = A„ 



a = 1,2. 



(69) 



When the two bands are anisotropic, wc can look for 
solutions of the self-consistency system (68) which satis- 
fies at Hc2 the linear equation 



-(e^)ifcninfcA„ = A„ 



a = 1,2. 



(70) 



All components of the tensor (^^)ifc are to be determined 
from the self-consistency equations. One can consider 
Eq. (70) as an ansatz which should be substituted in the 
self-consistency relations. If one succeeds in finding such 
a (C^)jfe so that the latter are satisfied, the ansatz (70) is 
proven correct. 

Repeating the derivation of Section IV we obtain: 



(71) 



where ^ is the average coherence length related to the 
eigenvalues of = ^2^2/3 ^^^j ^2^ _ ^2^-4/3^ p^j._ 

ther, T] = p^v^/AS^^ and v±_ is given in Eq. (45) for H \\ c 



Zero determinant of the linear system (72) gives S,{t): 
^''{Intf -e^^t{XiiJ^+\22J2) + 5JiJ2 = G. (75) 

For a single band n2 = and An = 1 and one obtains 
Eq. (44). 

The order parameters A1.2 at i?c2, as solutions of the 
system (72), are determined only within an arbitrary fac- 
tor, whereas their ratio is fixed: A1/A2 = —012/011- 
This means that Ai and A2 at Hc2 have the same coor- 
dinate dependencies and, in particular, that they have 
coinciding zeros (this, of course, follows already from 
Eq. (70)). The gaps ratio is, in general, temperature de- 
pendent (for brevity, we use the term "gap" instead of 
"order parameter" although the latter is more accurate). 

Introducing a dimensionlcss field he according to 
Eq. (49) one rewrites Eq. (75) as an equation for hc{t): 



(Int)^ - 2/ic(niAiiJi -|- n2\22'^2) Int 

+ 4/i^(niAii + n2A22 - l)XiX2 = 

dsslntanh(st)(^ce " 



Jo 



(76) 
(77) 



where is given by Eq. (49) for the corresponding band, 
and we took account of Eqs. (66) and (67). As in the one- 
band case, this equation does not contain the anisotropy 
parameter 7 and can be solved numerically for hc{t). 
Equations of a structure similar to (76) have been em- 
ployed in studies of Hc2 in two-band superconductors.*'^ 
Given hc{t), one finds the upper critical field: 



-"c2,c - ^2y2 "'c\t) , 



(78) 



where vo is expressed in terms of the Fermi energy and 

the total density of states in Eq. (47). 

Writing the self-consistency condition for if _L c, we 
obtain Eq. (76) in which, however, hc{t) is now known 
and /ica in integrals (77) is replaced with fih.ai"') given 
by Eq. (51) for each band. Solving this numerically, one 
obtains "/{t). 

The case of T ^ Tc is treated as was done for a one- 
band situation: 



7C(3) 
16 



(79) 



(take Eq. (53) , set J7 = 1 for the s-wave and add the band 
index). For -ff || c one takes ^ = fic whereas fi = nt, for 
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H J- c. The same argument which led to Eq. (56) gives 
for two bands at low temperatures: 



2a (m) 



2hc 



(ln(/tcM))c 



(80) 



One can make progress analytically in looking for Hc2 
near Tc and for T — > 0. This calculation is useful for 
checking the numerical routine; for the sake of brevity 
we do not provide these somewhat cumbersome results. 



A. Ratio A2/A1 at i?c2 

It follows from the system (72): 

Ai _ 022 _ Int - 2n2A22/ic22 
A2 a2i 2niAi2/ic2^i 

where 



Ai 



nin2AiiA22 — niXu — n2\22 + 1 

nin2 



(81) 



(82) 



We stress again that the coordinate independent ratio 
A2(r)/Ai(r) makes sense only for order parameters hav- 
ing the same coordinate dependence (in particular, the 
same zeros and the same phases). 

As T ^ 0, one can keep only the leading logarithmic 
term in of Eq. (80) to obtain:"*^ 



Ai 
A2 



1 — ri2A: 



22 



n2Ai 



(83) 



We point out that the zero-T gaps ratio does not depend 
either on the Fermi surfaces involved or on the value of 
ifc2(0). Also, this ratio at t = is the same for both field 
orientations; this is not the case for t^Q. 



VI. TWO BANDS WITH GAPS OF DIFFERENT 
SYMMETRIES 

Other than s-wave order parameters emerge if the cou- 
pling F(fe, k') responsible for superconductivity is not a 
constant (or a 2x2 matrix of k independent constants). 
The formally simplest way to consider different from s- 
wave order parameters without going to details of micro- 
scopic interactions is to use a "separable" potential: 



v^p{k,k') = v^l^n^{k)nf,{k') 



(84) 



where V^^ is a fc independent matrix, and look for the 
order parameters in the form 



A„(r,T,fc) = \E',(r,T)f]„(fc) 



(85) 



with the normalization (^^^)^ = 1 for both bands. One 
can see that this leads to the self-consistency equation 



-*c<lnt = 27rr^ 



Y,npKp{^pfp)^ I , (86) 





where 



(87) 



see also Appendix B. The same algebra as in Section V 
results in Eq. (76) for hcit) with 



^0^ = J rfss In tanh(st){^f2Vce"'''=''°^' 



(88) 



As before, one calculates the anisotropy 7(i) replacing /Xc 
with fibi-y)- 



VII. ELLIPSOID OF ROTATION 

The Fermi surface as an ellipsoid of rotation is an inter- 
esting example! on its own right and as a model system for 
calculating Hc2 in uniaxial materials with closed Fermi 
surfaces. Since Hc2 is weakly sensitive to fine details of 
Fermi surfaces, calculations done for ellipsoids might be 
relevant for realistic shapes as well. 

Similarly, open Fermi surfaces (extending to bound- 
aries of the Brilouin zone) in uniaxial materials can 
be studied qualitatively by considering rotational hyper- 
boloids. The formal treatment of these shapes is similar 
to that of ellipsoids. This work is still in progress, we 
show some of it in Appendix D. 

Consider an uniaxial superconductor with the elec- 
tronic spectrum 



Eik) 




(89) 



so that the Fermi surface is an ellipsoid of rotation with 
z being the symmetry axis. 
In spherical coordinates (fc, 0, (p) we have 



E{k) = sin^ e + ^ cos^ 0) = , 

2mab V J 2mab 



so that 



2 _ 2mabEF 



(90) 



(91) 



The Fermi velocity is v{k) = ^kE{k), with the deriva- 
tives taken at fe = feir: 

Vab sin 9 cos 4> Vab sin 9 sin (p 



vm 



Vab cos 9 rUab j^Ep 

The value of the Fermi velocity, v = {v'^ + Vy + vlY^"^, is 
given by 



V = Vab\ 



I SVC? 9 + 8^ cos^ 9 
s\v? 9 + e cos^ 9 



Vab\ 



m 



(93) 
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The density of states N{0) is defined as an integral 
over the Fermi surface: 



7V(0) 



{2TThfv 



(94) 



The average over the Fermi surface can be written as 
average over the solid angle dVL = sin 9 dO d4>: 



N{Q) 



dn 



2TT^h^ J 47rVrWW)' 
The Fermi surface average of a function A{6, 
1 f dflA{e,(j)) 



(A(0,0)> 



D = 



dn _ F(cos-i ye, 1 - e) 



(95) 

(96) 
(97) 



where F is an Incomplete Elliptic Integral of the first 
kind. If the function A depends only on the polar angle 
9, one can employ the variable u = cos 9: 



1 



du A{u) 



D{e) 7o ^T{u,e)T^{u,e) 
^=l + (£-l)u^ Fi = 1 + (e^ - 1)^2 



(98) 
(99) 



It is useful for the following to have a relation between 
Vab and uq of Eq. (47) for a one-band situation: 



vlb = D{e) vl 



(100) 



H 



One obtains /Zc using Eqs. (47), (49) and (D8): 



V{9,e) 



(101) 



Hence, we can solve Eq. (48) for hc{t) for a spheroid with 
a fixed mab/nic = e. 

Examples of this calculation for the s-wave order pa- 
rameter, $7 = 1, are shown in the upper panel of Fig. 1 
for a prolate ellipsoid with e — 0.2 and for an oblate one 
with e = 5 (the latter corresponds to the ratio of spheroid 
semi- axes -y/S); the numerical procedure is outlined in 
Appendix C. The plotted hc(t) is Hc2,c normalized on 
(j>o/{h'^VQ/27TT^) with vo given in Eq. (47) in terms of the 
Fermi energy and the total density of states. For com- 
parison, the same results are shown in the lower panel of 
Fig. 1 in the traditional HW normalization 



K = 



Hc2AT) _ hcjt) 



TcK2ATc) 



(102) 



It is seen, therefore, that for one-band s-wave materials, 
although the actual values of /ic(0) vary, the curves of 
Hc2,c{T) have qualitatively similar shapes for different 
Fermi surfaces. 




0.4 0.6 
t = T/T 

FIG. 1. (Color online) The upper panel: reduced upper criti- 
cal fields for a prolate e = 0.2 and oblate e = 5 spheroids and 
s-wave order parameter. hc(t) is calculated solving Eq. (48), 
(49), (52). The HW result for e = 1 is shown for comparison. 
The lower panel: the same result plotted using the HW nor- 
malization (102) to show that in this representation he only 
weakly depends on the Fermi surface shape. 



B. fit) 

To solve Eq. (50) for i? || a, we need 



^J■b 



=I?2/3(£) 



sin^ 9 cos^ ip -f f'^e^ cos^ i 



r(0,e) 



(103) 



The anisotropy parameter 7 is calculated with the help 
of Eq. (50), (51) and shown in Fig. 2. It is worth observ- 
ing that for the s-wave case, 7 depends on the Fermi 
surface shape but is temperature independent. 

One can see that 7 = 1 for a Fermi sphere with e = 1, 
as is should be. One can show that 7(e) behave approxi- 
mately as Xj \fe. In particular, we observe that for oblate 
Fermi spheroids, 7 < 1, i.e., Hc2,ab < Hc2.c- 



C. d-wave on a one-band ellipsoid 

For this case i7 = fio cos 2(p and one can verify that the 
condition (fi^) = 1 yields il^ = 2, the same value for any 
Fermi spheroid. Eq. (48) then can be solved numerically 
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e = 0.2 


1 1 1 1 ' 1 1 1 
,d-wave 




s-wave 
~ ' 


- 


- 

.d-wave 


■ e= 1 / 


' s-wave 


/ 




d-wave 




/^s-wave 


e = 5 

1 . 1 . 1 . 1 . 1 



0.0 0.2 0.4 0.6 0.8 1.0 
t 

FIG. 2. (Color online) Anisotropy parameter 7 = habile = 
Hc2,ab/ Hc2,c Calculated solving Eq. (50), (51), and (77) for s- 
and d-wave order parameters. 

with the results shown in Fig. 3 for values of e given in the 
caption. The anisotropy parameter for d-wave is shown 
in Fig. 2; unlike s-wave, it decreases on warming. 

It is worth observing that for a fixed e, hc{t) for s- 
and d-wave order parameters are the same. This is 
because the Fermi surface average in Eq. (48) involves 
{1/2) sniOdO for the s-wave whereas for the d-wave 

we have {I/Att) sin 6 dO J^^ 2 cos^ 2(j) d(j) which has the 
same value as for the s-case. 




t 

FIG. 3. (Color online) Reduced fields h{t) for two spheroids, 
e = 0.2,5. For H || c, s- and d-curves coincide, as explained 
in the text. 



D. Order parameter modulated along 

The gap function suggested by the ARPES data for 
Bao.6Ko.4Fe2As2 has a general form:"* 



so that 



A = Ao(l + ?7 cos k^a) . 



(104) 



This order parameter varies along the Fermi surface with 
changing fc^; it does not have zeros if \r]\ < 1. Depending 
on the sign of r/, it has maximum or minimum at the 
"equator" = 0. 

To apply this dependence for Fermi spheroids, we 
write: 



,2 ,2 2 a '^mabEp cos^ 9 



(105) 



where kp is taken from Eq. (D4). Since m|| = 2Ef/v'^^ 
and vn = vqD^/'^^ we obtain 



k,a = 



2Epacos6 



We now choose the length scale 

hvo 
''=2E^' 



(106) 



(107) 



cos 9 



i?i/3(e)yrM) 



(108) 



Note that in the isotropic case, the length a = 
00/(377^)^/"^ with gq being the interparticle spacing (the 
unit cell size). 

To adopt the order parameter (104) for our formalism 
we define 



(1 + ricos{kza)Y 
((1 + ricos{kza)Y) 



(109) 



as to satisfy the normalization (51^) = 1 (the average in 
the denominator is calculated according to Eq. (DIO)). 

Numerical evaluation of hc{t) for parameters given in 
the caption of Fig. 4 result in a curve qualitatively similar 
to that of HW (the upper panel), however, the anisotropy 
parameter decreases on warming for 77 > as shown in 
the middle panel. It is quite remarkable that changing 
the sign of 77, i.e., going from order parameters with a 
maximum at the Fermi surface "equator" fc^ = (77 > 0) 
to ones with a minimum {rj < 0), not only changes the 
temperature dependence of 7 to the opposite, but affects 
its absolute values as well (for positive 77, the anisotropy 
parameter is noticeably larger than for ry < for other 
parameters kept the same). 
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8= 0.1 




5.5 



5.0 
2.2 



ri = +0.5 



H \ < \ < \ < \ < h 



ri = -0.5 



0.0 0.2 0.4 0.6 0.8 1.0 



t 



FIG. 4. (Color online) The reduced field hc{t) is shown in 
the upper panel for a spheroid with e = 0.1 and the order 
parameter (104) with a — 1 and rj — ±0.5. Two lower panels 
show the anisotropy parameter j{t). 



One can readily evaluate how the anisotropy of the 
London penetration depth A, 7a = Ac/Aqj,, changes with 
temperature for -q — —0.5 (for which the Hc2 anisotropy 
is shown in the lowest panel of Fig. 4). To this end we 
note that the GL theory requires the same values of these 
two anisotropics at Tc, so that 7a (Tc) — j{Tc) ~ 2.17 ac- 
cording to Fig. 4. As T ^- 0, 7^(0) {vD / {vfi (in 
the clean limit, the order parameter does not enter the 
anisotropy of the London depth). The calculation 
of these averages is straightforward for a spheroid with 
£ = 0.1: 7a(0) « 3.38. Thus, the A-anisotropy de- 
creases on warming, unlike iJc2-anisotropy. This qual- 
itative behavior of the A-anisotropy is, in fact, seen in 
experiments on Ba(Fei_a;Co2,)2As2, (Bai_j;Ka;)Fe2As2, 
and NdFeAs(0i_^F^).4i 

In Fig. 5 we demonstrate that the effect of 7 increasing 
on warming remains if the Fermi surface changes from 
a prolate spheroid with e = 0.1 to a sphere e = 1 and 
to oblate spheroid with e = 5. In particular, these fea- 
tures challenge the common belief that temperature de- 
pendence of the anisotropy parameter is always related 
to a multi-band situations. 

We note in concluding this Section that other possible 
anisotropic order parameters can be treated within our 
scheme in a similar manner. 
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FIG. 5. (Color online) The anisotropy parameter 7(f) for 
three different Fermi surfaces e — 0.1, 1,5 and the order pa- 
rameter of the form (104) with 77 — —0.5. 



VIII. TWO-BAND RESULTS 

To apply the theory developed for two-band materi- 
als one first should map actual band structure upon two 
spheroids, the procedure we demonstrate in some detail 
on the well-studied MgB2. When calculating parameters 
fic (and fib) needed in this mapping for each band, one 
should bear in mind that in the two-band situation we 
have: 



see Appendix C. 



Djeg) 
n„ 



2/3 



sm a 

r(0,£a) 



(110) 



A. MgBa 

We take this example to demonstrate that our proce- 
dure yields Hc2,c{T) and the anisotropy 7(T) in agree- 
ment with existing data (see, e.g. Refs. 12 and 42) and 
with calculations of Refs. 2, 3, and 9. We stress that our 
calculations of Hc2 are done with the same set of coupling 
parameters as those used for the zero-field properties of 
this material as described in Ref. 39. The four Fermi 
sheets of MgB2 can be grouped in two effective bands 
with nearly constant zero-field gaps for each group. The 
two effective bands are mapped here upon two ellipsoids.^ 
We describe this procedure in some detail. 

We take the following data from the band structure 
calculations:''^ the relative densities of states of our model 
are ni « 0.42 and 712 « 0.58 for a and tt bands, 
respectively.^'^''** The band calculations*'* provide also 



the averages over separate Fermi sheets: 



23, 



{vl)i = 0.5, and (<)2 = 33.2, {vl)2 = 42.2 x 10" cm^/s^. 
To map this system onto two Fermi ellipsoids, we note 
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write: 
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FIG. 6. (Color online) The upper panel: reduced upper criti- 
cal fields of clean MgB2. hc{t) is calculated solving Eq. (76), 
(77) for two spheroids with ellipticity parameters ei — 0.02867 
(a strongly prolate spheroid for nearly 2D a band) and 
£2 ~ 1.273 (for a 3D vr band) evaluated with the help of 
two-bands Fermi velocities from the band structure calcula- 
tions of Ref. 44. The reduced couplings A^^ are evaluated on 
the basis of microscopic calculations of Ref. 45 as described in 
the text. The middle panel: the same for the fields in com- 
mon units; -ffc2,c(0) ^ 2.8 T and Hc2,ab{0) ~ 16 T. The lower- 
panel: the anisotropy parameter 7 — hat /he = Hc2,ab/Hc2,c 
calculated solving Eq. (76), (77) with /ic replaced by fJ^tif)- 



that the averages over spheroids are given by 



du{l 



2D{e) Jo r3/2(e)r;/'(£) ' 



9 9 



du \ 



D{e) Jo r3/2(e)ry^(e) 



(111) 
(112) 



where u = cos 9. The integrals here can be expressed 
in terms of Elliptic Integrals, alternatively they can be 
evaluated numerically. Forming the ratio we 
obtain an equation which can be solved for e. This gives 
El = 0.02867 and £2 — 1-273. For a given e and, e.g., 



(v^), we obtain t;"^ for two ellipsoids: 



3b,l 



2El 



D( 



(113) 



where Eq. (100) has been used; this gives Uq = 3.867 x 
lO^"' (cm/s)^, the constant used in the field normaliza- 
tion. 

To obtain normalized coupling constants Xap = 
Vajs/Vo we use the effective values (calculated including 
Coulomb repulsion) which in our notation read: iViVii — 
0.807, N2V22 = 0.276, N1V21 = 0.118, ^¥2^12 = O.OSO.-^^ 
Using Eq. (64) one evaluates l/VnN{0) = 1.211, and ob- 
tains the normalized coupling 



All = 



Vn ^ NiVn 1 
Vo N{0)Vo n. 



2.328. 



(114) 



Similarly, we find A22 = 0.5765, A12 = A21 = 0.340. With 
these input parameters we have solved Eqs. (76), (77) for 
hc{t). The result is shown in the upper panel of Fig. 6. 

Given hc(t), we rewrite the same equations where fic is 
replaced with Ubij) and solve them for 7(t). The latter 
is shown in the lower panel of Fig. 6. 

Hence, the calculation gives hc{0) ~ 0.292. We now 
use the relation (78) between the dimensionless he and 
physical Hc2,c, where Tc — 39.5 K and vq has been esti- 
mated above, to obtain i?c2,c(0) « 2.8 T, the value close 
to the observed. ^^'^^ The general behavior of Hc2,c{T) is 
close to that of HW, the fact confirmed by experiment. 
It should be mentioned that our result is close to the cal- 
culations of Miranovic et al,^ Dahm and Schopohl,^ and 
Palistrant.^ We, however, do not reproduce the calcu- 
lated Hc2,c{T) with changing curvature and a substantial 
upturn at low T's of Ref. 8. 




0.6019 X 10" and w^^^ ^ = 1-436 x 10^^ (cm/s)^ Next, we 



FIG. 7. (Color online) The ratio A1/A2 of two order param- 
eters of MgB2 at Hc2,c and at zero field calculated with the 
same coupling constants. 

The general formula for the gaps ratio at Hc2 is given 
in Eq. (81). Figure 7 shows this ratio as a function of tem- 
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perature at -ffc2,c- For comparison we show the gaps ratio 
in zero field calculated with the help of the same coupling 
constants. It is instructive to note that the gaps ratio at 
Hc2 exceeds the zero-field value at all temperatures, that 
can be interpreted as a stronger suppression of the small 
gap by the magnetic field than that of the large one at the 
leading a band. At first sight, one would expect the two 
ratios to coincide as T Tc, which our results clearly do 
not show. Such an expectation, however, would not be 
justified: even when Hc2 — > 0, the superconductor in the 
mixed state differs from the uniform state by an extra 
magnetic field suppression of the order parameter. 

One often finds in literature the statement that a small 
gap in MgB2 is substantially or even completely sup- 
pressed by a large enough field. This would correspond 
to a substantial increase of A1/A2 or even divergence of 
this ratio at some field under Hc2- Our result, however, 
shows that even at the gaps ratio is finite and of the 
same order at all T's. We conclude that the full sup- 
pression of the small gap does not happen at any field 
H < Hc2- On the other hand, assuming (as was done 
in Ref. 2) that the gap ratio at Hc2 is the same as that 
calculated in zero field, is also incorrect. 
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FIG. 8. (Color online) The field hc{t) (the upper panel) and 
the anisotropy parameter 7(t) (the lower panel) for the order 
parameter of the form A = Ao{l+ri cos k^a) . The parameters 
chosen are: £1 = £2 = 0.1, ni — 0.5, a = 1, for three sets of 
ry's given in the figure. 



Near Tc we use I's of Eq. (53) and he from Eq. (115) 
to obtain: 



< |Ai 



Ai 

A2 



1 



111X12 V (f^^/ic)l 



ln2{fl'^IJ-c)2 



ni{rfi^. 



(119) 



This case is close to theoretical models of pnictides 
in which the interband coupling is assumed dominant. 
We consider here a limiting situation An = A22 = to 
simplify the algebra. Indeed, for the two field directions 
we have: 



{\xitf ~ AhlIl{^i,)I2{^ic) 

(Intf -4hlli{fib)Mt^b) 



0, 
0. 



(115) 
(116) 



The first equation here can be solved for hc{t). Since 
/ib depends on 7, the second gives an equation for j{t). 
The latter can also be written in a different form if one 
subtracts Eq. (115) from (116): 



Il{^J-c)^2ifJ-c) ^Il{^J.b)l2i^J-b) 



(117) 



Fig. 8 shows numerical solution for j(t) for a few rep- 
resentative parameter sets. It is worth noting the par- 
ticularly informative feature of Fig. 8: it shows that the 
anisotropy 7(i) is not necessarily a monotonic function 
of temperature. We note that numerical calculations of 
7(t) showing an extremum at < t < 1 are quite robust. 



1. A1/A2 

The general result for the gaps ratio (81) gives for 
All = A22 = 0: 



Ai _ \nt 
A2 2niAi2/ic2^i 



(118) 



since for An — X22 — 0, the normalized A^2 = l/nin2. 
As r — > 0, we can keep only the leading term ^ ln< in I 
of Eq. (56): 



Ai 
A2 



T=o V ni 



(120) 



It is instructive to note that the last relation in the 
form TiiA^ = n^Aj at T = suggests the equipartition 
of condensation energy between the bands (provided the 
only non-zero coupling is A12). 



2. Two bands with A(fej) 

This example is of interest because it may have impli- 
cations for understanding the behavior of Hc2{T) and, 
in particular, its anisotropy in Fe-based materials. In 
Fig. 9 we show hc{t) and ^{t) for two nearly cylindri- 
cal bands with order parameters modulated along kz ac- 
cording to Eq. (104). Modulations are characterized by 

= ??2 = —0.2; other parameters are given in the cap- 
tion. The main feature worth paying attention to is the 
anisotropy ^{t) which is increasing on warming for neg- 
ative ?7's. 

It is worth noticing that experimental anisotropy of 
Fe-based materials behaves qualitatively similar to that 
shown by the lowest curve in the lower panel of Fig. 9.^^ 
We, however, do not have enough information to fix the 
necessary parameters for realistic calculations (one needs 
partial densities of states, Fermi surfaces characterized 
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FIG. 9. (Color online) The field hc{t) (the upper panel) and 
the anisotropy parameter 7(t) (the lower panel) for the order 
parameter of the form A = Ao(H-?? cos k^a) . The parameters 
chosen are: £i = £2 = 0.1, m = 0.5, a = 1, for three sets of 
jy's given in the figure. 

separately for relevant bands for evaluating the geomet- 
ric parameters e's, and the order parameters A(fci?)). 
Hence, we take our results as having a generic qualita- 
tive value. In particular, among various possibilities we 
have considered, only the order parameters of the form 
A = Ao(l+77 cos kzu) combined with dominant role of the 
inter-band coupling yield the j(t) increasing on warming 
similar to what is seen experimentally. 



IX. SUMMARY AND CONCLUSIONS 

The upper critical field i?c2 and its anisotropy are 
among the easiest properties to examin when a new su- 
perconductor is discovered. This work provides a rela- 
tively straightforward scheme for evaluating the orbital 
Hc2{T) and its anisotropy 7(T) for single and two band 
uniaxial materials. We reproduce here the main points 
of our approach. 

The input parameters for two band materials are (i) 
the coupling matrix (or the normalized couplings 
Aq./3 and Tc), (ii) the symmetries of the order parameter 
on two bands given as Vtp{9^(l)), f3 is the band index, 
with the normalization (3), and (iii) the characteristics 
of electron systems (Fermi surfaces, averages of squared 
Fermi velocities, DOS'). 

The equation to solve for the reduced upper critical 



field hc{t) parallel to the c crystal axis of a two-band 
clean uniaxial material reads: 

(Inf)^ - 2hc{ni\iili + n2A222;'2) Int 

+ 4ft2(„,An+n2A22- 1)112:2 = 0, (121) 

Ip= j dsslntanh(s<){|0Vce"^'^"''''=^^ , (122) 

where is given by Eq. (49) for the corresponding band, 
and rj's describe the order parameter symmetries. After 
hc{t) is found, one solves Eq. (121), in which is re- 
placed by ^6(7) of Eq. (51), for ^{t). The one-band case 
is obtained by setting ^2 = 0,ni = 1 and An = 1. The 
case of two s-wave bands corresponds to ili = fi2 = 1. 

Because iJc2 is determined by equations containing 
only integrals over the Fermi surface, it is insensitive to 
fine details of the Fermi surface shape. Therefore, one 
can replace actual Fermi surfaces with ellipsoids (or with 
spheroids, for uniaxial materials). Given the averages of 
the squared Fermi velocities over each band, one estab- 
lishes the geometry of corresponding rotational ellipsoids 
(the squared ratio of semi-axes, e). This procedure is de- 
scribed in Section VIII. A on the well-studied example of 
MgB2. 

One numerically solves Eqs. (121), (122) by employing 
any of available packages (such as Mathematica) able to 
find roots of nonlinear equations and to do multiple inte- 
grals. The scheme can also be applied to the case of two 
bands with order parameters of different symmetries. 

By design, our method is applicable for clean mate- 
rials with a moderate -ffc2(0); paramagnetic limiting ef- 
fects are out of the scope of this work.® The method 
differs from those previously employed by not involv- 
ing explicit coordinate dependent A1.2 and minimiza- 
tion relative to the vortex lattice structures in calculating 
Hc2{T)P''^'^^'^^ The main feature of the two-band deriva- 
tion is that the linearized GL equation (70) is assumed to 
hold at Hc2{T) at any T, the ansatz proven correct by sat- 
isfying the self-consistency equation of the theory. The 
method is tested on the well-studied example of MgB2 
where it shows a satisfactory agreement with data and 
with other calculations. 

Our main results are as follows: 

1. We find that in clean one-band s-wave materials, 
the T dependence of the anisotropy 7 cannot be caused 
by the Fermi surface anisotropy (however, the param- 
agnetic limit, which is out of the scope of this paper, 
may suppress 7 at low T's and cause 7(T) to increase on 
warming) . 

2. For other than s-wave symmetry, 7 depends on tem- 
perature even for one-band materials. This dependence 
is pronounced for open Fermi surfaces as well as for or- 
der parameters depending on kz- Thus, the common be- 
lief that the temperature dependence of the anisotropy 
parameter is always related to multi-band situations is 
incorrect. 

3. Our scheme of calculating 11^2 for two-band ma- 
terials does not utilize any assumptions about the field 
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and temperature dependences of the order parameters 
Aa in two bands. ^ In fact, the gaps ratio is calculated 
self-consistenly and in general turns out temperature de- 
pendent. Although both Ai and A2 turn zero at Hc2{T), 
their ratio is finite and in the examples examined is larger 
than at zero field. 

4. The case of exclusively inter-band coupling is dis- 
cussed, that might be relevant while interpreting data on 
Hc2 and its anisotropy in Fe-based compounds. 

5. For order parameters of the form A = Ao(l -t- 
?7cos kzo) (one of the candidates suggested for pnictides), 
the anisotropy parameter j(t) depends on the sign of 77 
(or ?7's for two-bands). In particular, "/{T) increases on 
warming in a nearly linear fashion (as for pnictides) for 
both rj's negative. 
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Appendix A: Different form of the one-band 
equation for Hc2 



Appendix B: Tc as a function of Aa/3 

This question had been discussed, e.g., in Ref. 39). 
Since our notation of normalized X^p differs from that in 
the literature, we provide here corresponding relations. 
The s-wave self consistency equation for H = 0, 



27rTV7V^K^/^(a;) 



gives near where = A^/?ia;: 



2'kT 



where Vq is to be defined. We choose it so that 

1 ■^S^ 2'kT ^ 2fiujD 

N{0)Vo ~ ^ luv " ^ 1.76Tc ' 



(Bl) 



(B2) 



(B3) 



X^ju J^AiAi^juA^, zero determinant of which 



or, which is the same, l.l&T^ = 2hwDe~^/^^^^^° . We 
then obtain a linear and homogeneous system of equa- 
tions Ajy 
gives Eq. (66). 

For other than s-wave order parameters on two bands, 
we take the coupling potential in the form (84) and the 
order parameters as in Eq. (85) , We then obtain the self- 
consistency equation = '2''^TY,^ ,^N^v!Pj{Q.^f^{ijj)). 

We now denote X^^ = vJ^V^o, Eq. (87), and recall that 
ff, = ^lf,'^f,/huj near T,.. This gives *y = I]^ ^^^Aj^^^"^, 
i.e. the same system of equations as above for the s-wave 
case and the same zero-determinant condition (66) albeit 
with renormalized couplings (87). 



Both sides of Eq. (48) diverge logarithmically when t 
0, so that these divergences, in fact, cancel out. However, 
in numerical work this cancellation may not always be 
exact which may cause the numerical solutions to become 
unreliable in this limit. Here, we provide an alternative 
form of this equation free of this shortcoming. 

To this end, consider an identity: 

\nt = (n'^2hciJ. dssln{st)e- ''''''' 
C + {\n{n''hafi)) 



+ 



(Al) 



which is verified by direct integration. We now combine 
it with Eq. (48): 



C + (1^2 ln(/ie/x)) 



Jo 



ds s In 



(A2) 



As t 0, the integral on the RHS goes to zero, and we 
immediately obtain the result (57) for hc{Q). 



Appendix C: iXc for the two-band case 

By definition, ji^p oc v'^^ ^/vq with vq given in Eq. (47). 
For Vab,p we have two relations: one with the Fermi en- 
ergy, mabj/aVai, = 2Ef, and the other in terms of the 
band DOS, Eq. (95): 



27r2^3 



(CI) 



One excludes mab,p from these two relations to obtain: 
Hence, we have: 



Mc,/3 



^ab,l3 



Die, ) 



1 2/3 



(C3) 



a clear generalization of Eq. (100) for the one-band 
spheroid. This gives Eq. (110) used in the text. 
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Appendix D: Open Fermi Surface 

The theory employed above is designed to model Fermi 
surfaces closed within the first Brillouin zone. Here we 
consider an example of the Fermi surface which crosses 
the zone boundary, i.e. it is open. Perhaps, the simplest 
shape to consider is a rotational hyperboloid which is a 
property of the carriers energy of the form: 



E{k) = 




(Dl) 



A schematic picture is shown in Fig. 10. In spherical co- 




FIG. 10. (Color online) The cross-section fc^ = of an open 
Fermi hyperboloid. kz,BZ is the zone boundary in the c di- 
rection. Tilted straight lines correspond to -E = 0, whereas 
the hyperbola shows the Fermi surface. 



ordinates (fc, 6 
E{k) ^ ^ 

^2 



2mab 



b) we have 
sm 



cos 



ecos 



rric J 2mab 

e = mab/mc . 



^2{0), 



(D2) 



It is seen from the figure that in the first quadrant of the 
plane kx,kz, the Fermi surface is situated aX 9 > 9m > 
tan~^ y/e, i.e., everywhere at the Fermi surface F2 > 0. 

The angle 9^ corresponds to the crossing of the Fermi 
hyperboloid with the zone boundary kz^sz = 27r/c where 
c is the unit cell size along the c direction: 



tan 9m — e 



mabC^EF 



kl{^/2)/kl 



= e + a. (D3) 
in most situations of 



The parameter a — n,pyn / z,; / n,^ 
interest is less than unity {kp{TT/2) is the radius of the 
hyperboloid neck). 

The Fermi momentum is given by 



2mabEF 
h^T2{9) 



(D4) 
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FIG. 11. (Color online) The upper panel: /ic(t) for Fermi hy- 
perboloid with e = a = 0.1. The lower panel: the anisotropy 
7(i) for the same parameters of the Fermi hyperboloid. 



The Fermi velocity is v{h) ~ V kE{k)kp: 



Vab sm a cos ( 



Vab sm a sm ( 



Vab COS 9 



Vab = \ 2EF/mii 



Futher, we have for v = {v^ + + ^zY^^'- 



(D5) 



V = Vab\ 



I sin^ 9 + 6^ cos^ ( 
sin^ 9 — e cos^ 6 



Vab\ 



^2{9) 



(D6) 



The density of states N{0) is defined by the integral 
of Eq. (94), which can be written as an integral over the 
solid angle dfl — sin 9 d9 dip: 



NiO) 



m-ab^ab 



dn 



2Tr^h3 J 47rVrW)rW) 



(D7) 



where the integration over 9 is extended from 9m to tt - 
9m- The Fermi surface average of a function A{9, (j)) is 



D = 



1 



dnA{9, 



Dj A7T^r2{9)r,{9) ' 
dn 

For A depending only on 9, one can employ u = cos 
1 /■"'" duA{u) 



(D8) 
(D9) 



{A) = 



(DIO) 



Die) Jo VF2(u,e)Fi(u,e) ' 
F2 = 1 - (e -t- 1)^2 , Fi = 1 -f (e^ - 1)^2 . (DU) 
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where the upper hmit is Um = cos 9m = l/\/l + e + a. 
In particular, one obtains: 



Vl + e 

where F is an Incomplete Elliptic Integral of the first 
kind. 

As for ellipsoids, the relation between Vab and vq de- 
fined in Eq. (47) for a one-band situation holds: 



vL = Dv? 



; 



(D13) 



however, with a different D. 

For H II c, the relevant electron orbits are circular, as 
for spheres and rotational ellipsoids, and we do not ex- 



pect qualitative deviations from the latter. Fig.?? shows 
hc{t) calculated with the help of Eqs. (48) and (49) for 
both s- and d-wave order parameters. We note that the 
anisotropy parameter shown in this figure is suppressed 
substantially on warming. 



Numerical evaluation of the anisotropy shows that 
7(t) decreases on warming. Fig. 11. Reduction of the 
anisotropy on warming is substantial for a single band 
open Fermi surface, an interesting observation given 
the common belief that the temperature dependence of 
anisotropy is a multi-band property. We also note that 
this reduction is stronger than that of open Fermi sur- 
faces. 
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